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Abstract: Precise measurements of SM particles properties at the LHC allows to look for 
heavy New Physics in the context of an Effective Field Theory (EFT). These searches, how¬ 
ever, often rely on kinematic regions where the validity of the EFT may be compromised. 
In this paper we propose to address this issue by comparing with benchmark models. The 
connection between models and their manifestations as EFTs at low energies allows us 
to quantify the breakdown of the EFT, and describe ways to combine different sources of 
constraints beyond Higgs physics. To illustrate these techniques, in this paper we propose 
a set of benchmark models based on extensions of the Higgs sector, namely the inclusion of 
a singlet, a dilaton and generic 2HDMs. We obtain the matching between these models and 
the EFT involving the Higgs, electroweak bosons and fermions. We then describe current 
and future indirect and direct constraints, consider the effect of correlations among the 
coefficients within models, and discuss the validity of the EFT. 
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1 Introduction 

The discovery of a Higgs boson [1, 2] has opened the Era of Higgs Physics. An intense effort 
is now devoted to measure the properties of this new particle h and determine whether 
it is the Standard Model (SM) Higgs boson [3-6]. So far, no evidence for new physics 
beyond the SM has been observed, which suggests the use of an Effective Field Theory 
(EFT) approach, that assumes possible new physics beyond the SM modifying the Higgs 
interactions to be heavy, with a typical scale A ~ TeV. The effects of new physics are then 
parametrized through higher-dimensional operators constructed from the SM fields, the 
leading operators appearing at dimension-six [7-12]. Some of these operators which affect 
Higgs properties can be measured in Higgs physics only, while others are also related to 
electroweak (EW) observables since the Higgs scalar excitation is always associated with 
the EW symmetry breaking order parameter v. The experimental data from LEP and 
Tevatron constrain the size of the Wilson coefficients of these operators, and more recently 
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the results from the LHC experiments ATLAS and CMS allow also to constrain the Wilson 
coefficients associated to the operators that impact Higgs physics [13-21], 

The LHC energy reach allows to go beyond the analysis of Higgs signal strengths as a 
probe of new physics beyond the SM, exploiting the information encoded in kinematical dis¬ 
tributions to boost the sensitivity to new physics beyond the SM in the Higgs sector. This 
has the power to strongly constrain the presence of new physics which lead to a significant 
enhancement of the signal in certain kinematical regions, typically at high energy. How¬ 
ever, while the bounds on the Wilson coefficients of these dimension-six operators which 
can be extracted from these measurements are indeed very strong even with the limited 
amount of present data, the validity of the EFT approach in the kinematical regions which 
dominate these constraints is questionable [22, 23] (see also [21]). This would render these 
constraints meaningless in the context of EFTs. It is therefore crucial to assess in detail the 
validity of this approach, and a way to do so is by quantitatively studying the breakdown 
of the EFT by comparing its predictions with those of possible UV completions. 

Among possible UV completions, extended scalar sectors provide an attractive arena 
for the use of a SM effective theory description. Extensions of the SM Higgs sector consti¬ 
tute a simple yet very well-motivated scenario beyond the SM, with important consequences 
not only for phenomenology but also for EW cosmology, baryogenesis and dark matter. It 
is plausible that the new scalar states, which we assume to be significantly heavier than the 
Higgs boson discovered by ATLAS and CMS, cannot be detected directly at LHC either 
due to them being very heavy or having vanishingly small couplings to SM gauge bosons 
if the light Higgs is SM-like. However, it may be possible to measure the effect of these 
new particles via EW precision observables (EWPO) or using the kinematical information 
in LHC measurements. An EFT of extended Higgs sectors provides a very useful tool to 
study these effects systematically under the assumption that the new states are signifi¬ 
cantly heavier than h. Moreover, since extensions of the SM Higgs sector provide a very 
simple UV completion to such an EFT, this makes it possible to compare the predictions 
and bounds derived from the EFT with those of its UV completion, therefore probing the 
range of validity of the EFT. This allows us to assess the reliability of the LHC constraints 
drawn from the high-energy kinematical regions, as well as those from EWPO, and also to 
analyze the LHC potential for indirectly probing extended Higgs sectors. 

The paper is organized as follows: In Section 2 we introduce the D = 6 effective 
operators which are relevant for the analysis. In Section 3 we present the matching of 
these operators with UV models corresponding to extensions of the SM Higgs sector with 
an extra scalar singlet, with an extra scalar doublet (a Two-Higgs-Doublet-Model) and 
with a radion/dilaton, and discuss their experimental constraints both from the UV model 
and from the EFT point of view. In Section 4 we discuss the validity of the EFT in 
the light of these scenarios. Finally, in Section 5 we summarize our results and discuss 
their implications. In the Appendices we describe the connection to the various possible 
D = 6 Lagrangian terms after EW symmetry breaking, parametrized as anomalous Higgs 
couplings, and give details on the 2HDM relevant to our analysis. 
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2 Effective Field Theory for the SM Higgs Field 


The Lagrangian for the Standard Model can be supplemented by higher-dimensional op¬ 
erators that parametrize the effects of new physics beyond the SM appearing at energies 
much larger than an effective scale which may identified with the vev of the Higgs field v. 
Considering only operators of dimension D < 6 and assuming baryon and lepton number 
conservation, the most general SU(3)c x SU(2)l x t/(l)y gauge invariant Lagrangian £es 
constructed out of the SM fields has been known for a long time [7-9], and may be mapped 
into various alternative bases of D < 6 independent SM effective operators [10-12] (see 
also the discussion in [24]). 

Each experimental measurement will in general constrain only a handful of effective 
operators of a certain basis. In combining the available constraints from analyses of various 
physical processes (via e.g. a global fit), some basis choices will be more appropriate than 
others 1 . As compared to those of [10, 12], the basis of operators Oj from [11] is well- 
suited for analyses of Higgs properties in combination with precision measurements of EW 
observables. In this basis, the SM effective Lagrangian £e& may be expressed as 

^Efi = £sm + ^ CiOi = £sm + £silh + £>g + £cp + £f 3 + £f 2 + £f 3 ■ (2.1) 

i 

We adopt here the decomposition in [25, 26] and normalize the Wilson coefficients c* ac¬ 
cordingly. The term £silh of (2.1) corresponds to a certain set of CP-conserving operators 
involving the Higgs doublet 4>, inspired by scenarios where the Higgs is part of a strongly 
interacting sector [11] 

£silh = - ^[fU] 3 

,n w 

+ 219 C 2 HW [Di*&T 2k D v §\ W* u + [Df*&D v §\ B, IV (2.2) 

'W'W 

+ 9 -^&<S>B^B^ + GT 

- '^&$Q L &u R + ^&$Q L $d R + ^-&$L L §l R . 

L v z v z v z J 

Here A stands for the Higgs quartic coupling, g ', g and g s are respectively the t/(l)y, 
SU(2)l and SU(3)c coupling constants and T 2k = (Jkt^ axe the generators of SU(2) in 
the fundamental representation (a k being the Pauli matrices). The Hermitian derivative 
operator D M is defined as 


& D ^ = & D tJ (L> - , 


(2.3) 


1 Depending on the purposes, bases different from the one we adopt here might be more adequate. 
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and our conventions for the gauge-covariant derivatives and field-strength tensors (following 
[26]) are 


= d,jB v - d u B p 

= d,Wt - d v Wi + g e. h k W* W> 

G% = d.Gl - d v Gl + g s f bc a Gl Gl 

= d^-ig T 2k W k <$> - l -g' B p <S> (2.4) 

D p W k v = d p d p W k - d p d u w k + g ej d p (Wj W*) 

+g ej W; ( d p Wl - dyWj L ) - g 2 W pi (wj W k - W* Wj) 

D p G“ u = d p d p Gl - d p d v Gl + g s f bc a d p (Gj G d ) 

+9s f bc a G b p ( d p Gl - d u G c p ) - g 2 * s G pb (g\ G a v - G b G“) , 


with Cjj'’ and f bc a being respectively the structure constants for SU(2) and SU( 3). The 
term Cq consists of operators not directly connected to Higgs physics, but that affect the 
gauge sector through modifications of the gauge boson self-energies and self-interactions, 


Cg = 9 ^e ijk Wi v W^ 0 W^ k + g -^fabcG a U yG^ p G^ c + jjpW k D p W p v 


m 


w 


mt 


mt 


+ 

mi, 


d 0 B pv + 9 ^^D p G“ u D p G pu 


(2.5) 


[IV^p 


m 


w 


The term Cqp in (2.1) supplements £silh and Cq with a set of CP-violating operators 


C C P = 19 °!l W Di J ^T 2k D^W k + l ^ C ^D p &D u <$>B llv + <§>B pv B pv 


+ 


o =- -L. ZK, J -' ’ r uis I O 1 1 [IV I 9 ^ [±V A 


rnf 


mt 


"w ""w 

with the dual field strength tensors defined by 


mt 


( 2 . 6 ) 


5 1 

JJfj.i' — 2 & a wp<j 


i 


B pa , W* v = 2 e, lupa W^ k 


a _ _ s^ipcra 

^pv ~ 2 t [ivpcr^ r 


(2.7) 


Finally, there are further operators contained in £p 1 -operators involving two Higgs fields 
and a pair of quarks/leptons-, Cp 2 -operators involving one Higgs field, a gauge boson 
and a pair of quarks/leptons- and Cf 3 -four-fermion operators- (see e.g. [25, 26]). These 
effective operators are nevertheless not present at leading order in the extensions of the 
SM we consider in the present work - ’, and so we do not discuss them in the following. 
Furthermore, we consider CP conserving scenarios for the time being, leaving an analysis 
of Ccp for the future [27]. 

2 An exception are certain operators in Cf 3 which do get generated in Two-Higgs-Doublet-Model exten¬ 

sions of the SM. These are however proportional to Y a 2 (with Y„ a fermion Yukawa coupling), being then 

negligible for 1 st and 2 nd fermion generations, and so are essentially unconstrained. 
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After EW symmetry breaking, we can write the SM effective Lagrangian £es in the 
unitarity gauge and in the mass basis, with 



1 


{Wl^iWD 






(2.8) 


with s w (c w ) being the sine (cosine) of the weak mixing angle at tree level. In Appendix B, 
we detail the relation of the D = 6 terms presented here with the possible anomalous Higgs 
couplings after EW symmetry breaking. 


We are now ready to move onto relating the EFT to specific UV completions. In this 
work we have chosen to focus on the matching of the Wilson coefficients onto models with 
extended Higgs sectors, where the effects of New Physics are more apparent in deviations of 
the couplings of the Higgs and electroweak bosons. We describe the results of the matching, 
as well as the relation between these models and the EFT in the next section. 


3 Higgs EFT from Extended Higgs Sectors 

As discussed in the Introduction, an important aspect of the D < 6 SM effective theory from 
the previous section is its energy range of validity. This is a key issue if one is to reliably 
confront the predictions of the effective theory with experimental data, in particular those 
sensitive to energy scales E 3> v. In this sense, extended (non-minimal) scalar sectors 
provide a very simple renormalizable completion to the D < 6 SM effective theory, and 
allow for a quantitative assessment of the EFT’s energy range of validity. 

Moreover, extensions of the SM scalar sector provide an attractive arena for the use of 
the SM effective theory: they are a simple scenario beyond the SM, well-motivated from 
the point of view of EW cosmology and baryogenesis, and may also be regarded as part 
of a complete theory beyond the Standard Model at the TeV scale, such as Composite 
Higgs scenarios or low-energy Supersymmetry. Assuming that the new scalar states are 
significantly heavier than m/j, the effective theory of extended Higgs sectors provides a way 
to study these effects systematically, as has already been shown e.g. in [30]. In particular, 
it is plausible that the new scalar states from the extended Higgs sector are very hard to 
probe at LHC (they might be very heavy, or have no decay branching fractions to SM gauge 
bosons), but their effect might be possible to detect either in precision EW measurements 
or using differential information in LHC measurements. 

Below we construct the D < 6 SM effective theory for various extensions of the SM 
scalar sector: a singlet extension of the SM (the so-called “Higgs portal”), a Two Higgs 
Doublet Model and an extension of the SM by a dilaton/radion. In each case, we obtain 
the Wilson coefficients for the D = 6 effective operators by matching to the UV theory, 
and perform an analysis of their current experimental bounds and future prospects. The 
results of this construction for the different scenarios is brieffy summarized in Table 1. 

In the respective matching procedure we demand that the effective action for the UV 
theory and the D < 6 EFT agree after an expansion in the light degrees of freedoms over 
the new physics mass scale. This results in the matching of the one-light-particle irre¬ 
ducible (1LPI) Green’s functions in the full and effective theory - see [31] for a pedagogical 
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introduction for these type of calculations. In the following we will perform an off-shell 
matching, where we expand in external momenta and the EW mass-scale over the new 
physics mass scale and keep equation of motion vanishing operators in the calculation until 
the final projection. Accordingly all propagators of SM fields will be massless after the 
expansion and we can perform the calculation in the SU(2)l x U(1)y symmetric phase, 
quite analogous to the matching calculation performed in [32]. The light degrees of free¬ 
doms of the SU(2)l X U(l)y symmetric phase comprise the gauge singlet B, triplet W a 
and octet G a fields as well as the scalar doublet and the fermionic doublets and singlets. 
In practice we compute 1PI Green’s functions with up to 6 Higgs and 3 gauge boson fields 
where the total number of fields does not exceed 8. The resulting expressions are obviously 
related via SU{2)l x U(1)y gauge invariance of the operators in 2.2, 2.5, 2.6 and 3.40, 
which provides a useful consistency check of our calculation. Moreover, we do cross-check 
with an explicit calculation in the broken phase, detailed in Appendix B. 



CH 

C6 

CT 

c w 

C B 

C-HW 

Chb 

C 3 W 

C-y 

C 9 

Higgs Portal (G) 

L 

L 

X 

X 

X 

X 

X 

X 

X 

X 

Higgs Portal (Spontaneous (fP) 

T 

L 

RG 

RG 

RG 

X 

X 

X 

X 

X 

Higgs Portal (Explicit (jP) 

T 

T 

RG 

RG 

RG 

X 

X 

X 

X 

X 

2HDM Benchmark A (cg_ Q = 0) 

L 

L 

L 

L 

L 

L 

L 

L 

L 

X 

2HDM Benchmark B (c^_ Q ^ 0) 

T 

T 

L 

L 

L 

L 

L 

L 

L 

X 

Radion/Dilaton 

T 

T 

RG 

T 

T 

T 

T 

L 

T 

T 


Table 1. Leading order at which the various Wilson coefficients for the D = 6 SM effective field 
theory are generated in each of the scenarios under consideration. In each case, the operator can be 
generated at Tree-Level (T) or 1-Loop (L). If some operators are generated at Tree-Level, this may 
lead to the generation of others via operator mixing under 1-loop Renormalization Group evolution 
(see e.g. [28, 29]), which we denote by RG. Operators which are generated at higher order in RG 
and EFT expansion are denoted with an X. 


3.1 The Singlet Higgs Portal: Doublet-Singlet Mixing 

The addition of a singlet (real or complex) scalar field is arguably the simplest possible 
extension of the SM. Despite its minimality, this extension of the SM can have important 
consequences for the stability of the EW vacuum at high energies [33, 34], and it could at the 
same time constitute a “Higgs portal” into a dark/hidden sector [35, 36]. It may also have 
important consequences for Cosmology, potentially accounting for the dark matter relic 
density [37-39] or yielding a first order EW phase transition in the early Universe [40-47] 
that could explain the matter-antimatter asymmetry in the Universe through baryogenesis. 
In addition, it may give rise to interesting collider phenomenology (see e.g. [48-51]). 
Altogether, the singlet scalar extension of the SM constitutes a well-motivated scenario, 
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and the interplay of different present and future experimental data to probe it has been 
widely studied (see [52-54] for up-to-date analyses). Let us then consider the SM scalar 
potential extended by a singlet scalar field s 


F($, s) = -nl |<L| 2 + A |4>| 4 - y s 2 + ^ s 4 + ^ |<L| 2 s 2 . (3.1) 

We assume initially that linear and cubic terms in s are absent from E(<3?,,s), which may 
be achieved by means of a discrete/continuous symmetry G in the hidden sector, and focus 
on the scenario in which the field s develops a vei? , s —> s + v s . This generates linear and 
cubic terms in s, which specific relations among these and the rest of parameters in the 
potential. After EW symmetry breaking, the scalar potential reads 

V(h, s ) — h 2 + -^s 2 + ml s hs + vX k 3 + v s X 3 s 3 + / s 2 _|_ h 2 s 

+ + T fc “ + (3-2) 

4 4 4 

with = 2Au 2 , m 2 = 2 A s w 2 and m 2 s = X m vv s , and where we have used the min¬ 
imization conditions = Xv 2 + ^v 2 and /j, 2 = X s v 2 + %W 2 in V(h,s ) to trade the 
mass parameters /r^, /r 2 for the vevs. The term m\ s hs in (3.2) induces doublet-singlet mix¬ 
ing, leading to two mass eigenstates the lighter of which (hi) we identify with the 
discovered 125 GeV Higgs particle. The mixing angle 8 and masses are given by 

s 2 = _ 4 "4s _ = _V_ 

4 m^ s + (m 2 - m\ + xj(m 2 - m|) 2 + 4mJ s j 4 y 2 + (l — x 2 + a/(1 - x 2 ) 2 + 4 y 2 ^ 

(3.3) 

m i,2 = \ (j n h + m ‘s T \J(m 2 s - m 2 h ) 2 + ^ (l + x 2 T V i 1 ~ x 2 ) 2 + A y 2 ^ 

(3.4) 

with sg = sin(0), x = m,h/m s ~ v/v s and y = m\ s /m 2 ~ v/v s . The limit m s S> v 
corresponds then to v s v, with i,t/<l and so 


Sg — 


2 2 / 2 2 \ 2 2 
m 1 — m s (x — y ) = m h — s 6 




m 2 2 ~ m 2 (1 + y 2 ) = m 2 (1 + sg) (3.5) 


where we have neglected terms of 0(x 4 , x 2 y 2 , y 4 ). From (3.2), the relevant scalar self¬ 
interactions h\ h2 and h\ read 


v(hiM) d 




/*? + 


m 2 s y 

2v 


+ <D(x,y) 


h\ h 2 


(3.6) 


3 The spontaneous breaking of a discrete or global continuous symmetry G would respectively lead to 

domain wall formation in the early Universe or the existence of massless Goldstone bosons, both features 
being undesirable in a realistic model. Possible solutions are to consider G to be a spontaneously broken 
gauge symmetry, or to allow for a small explicit breaking of the symmetry. We will disregard these issues 
in the following discussion. 
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Neglecting 0(y 4: ) corrections, the couplings of h \ : 2 to the W ^ and Z bosons read 


g m w (1 - y 2 /2) hi + (l - y 2 ) h\ - ( gm w y ) h 2 + ... 


w+w* 1 - 

H 1 


+ 


gm z 


9 

Z 
C W 

-.2 1^,2 




-(^y i 




(3.7) 


Noting that y = mj ls /rn 2 = m\ s /m\ + 0(y 3 ), we may integrate out the heavy, singlet¬ 
like state h 2 • At leading order, this generates an 0{y 2 ) contribution to h\ V^V 11 , which 
precisely cancels the 0{y 2 ) correction in (3.7). The Higgs-gauge interactions then read 


a 2 ' 

gm w (l - y 2 / 2 ) hi + -- h\ 


w+w 1 *- + 


gm z 

2cw 


(1 - y 2 / 2) hi + 7^rh\ 

O C- 


w 


Z M Z^, (3.8) 


leading to a g\y 2 1 2 mismatch (with gy = — for V = W, Z respectively) between 

,oX „„(3) 


gUhvv and 9 om VV • We turn now to the fermionic couplings, focusing on the Higgs-top quark 


2 m w 


interactions 


Y 


V2 


Yy 


— 7 = 1 —tt ) hi t L t R - -= h 2 tttR , 


V2 


(3.9) 


with Yt the top-quark Yukawa coupling. This affects the gluon fusion Higgs effective 
coupling for hi, which gets modified w.r.t the SM one 


- 9j f (i - y 2 /2) a^ai-h , 


(3.10) 


The coupling G a lw Ga'h 2 (relevant for di-Higgs production) also gets modified, both through 
an 0(y 2 ) correction due to the singlet-doublet mixing and directly via the presence of h 2 , 
which can mediate the process gg —>• hihi (for s /m 2 s <C 1, being s the partonic center of 
mass energy for the process, this contribution is however small). 

The above discussion may be directly mapped into an SU(2)l x U(1)y invariant ef¬ 
fective field theory for the SM. We consider (3.1) after s develops a vev 

V($, s ) = -Ah l^l 2 + a |H>| 4 + 7 -fs 2 + V S X S s 3 + ^s 4 + X m v s m 2 s + ^ |<f>| 2 s 2 , (3.11) 

with Ah = hi ~ (Am^ 2 )/2. We then integrate-out the field s, which yields the effective 
Lagrangian 


A 2 ^ 2 
2 mi 




3n$i 2 aj$r- 6 


AmAs't’s 


-3 


A;), v 2 


rrvi 


mi 


|$| 6 = A _ a /^|$| 2 a M |$| 2 . ( 3 . 12 ) 


From (2.2), this corresponds to ch = y 2 . The two contributions to c§ are shown in Figure 
1, and cancel by means of rn 2 s = 2X s vl (yielding cq = 0). Upon EW symmetry breaking, 
ch modifies the wave function of the Higgs h 


Cee 7 ) (1 + ch) — [d^h) 2 


(3.13) 


- 8 - 






















, $ 




$ 


$■ 


TS 

S v*.. S 


— -• ^m'Vs 

,jfr A s u s 






' $ 


$ 


A, 


\ / 
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Figure 1. The two Feynman diagrams leading to the effective operator Cq in the singlet extension 
of the SM. Their respective combinatorial factors are 6 (Left) and 3 (Right). 


which universally reduces all couplings of h: e.g. the coupling between h and the weak 
bosons simply read 


(3) 9v v 2 
9hVV = n -y - 9V v 

v 1 + v 



(3.14) 


matching as expected the result from (3.8). By comparing the EFT results with (3.8) we 
observe that the EFT appproach reproduces the effect of scalar mixing on interactions 
involving one Higgs scalar h, but fails to do so for the case of two scalars hh. This suggests 
that di-Higgs production (see e.g. [49-51, 55, 56]) could be used as a tool to discriminate 
between a pure linear realization of the SM Higgs doublet and the presence of mixing 
effects, although the prospects seem challenging and we will not discuss them here. 


The previous results for ch are left unchanged by the inclusion of explicit linear and 
cubic terms for s in (3.1). These can however generate a non-zero value for cq, since now 
the cancellation among the diagrams in Figure 1 is not exact. Let us illustrate this by 
adding a term (x m |4>| 2 s to (3.11). This term does not alter the minimization conditions, 
but contributes to the singlet-doublet mixing upon EW symmetry breaking, so that now 
the mixing is given by y = m^/m 2 = (A m u s + y, m )v/m 2 . Upon integrating-out s, this 
now yields 




(Am Vs T 9m)* 

2 mi 


^| 4 >| 2 9 m |$| 2 - ( 6 


(Am Vs T 9m)^^s'Us g'^m^mhs T 9m) 2 


m'i 


m ’ 


'%dn<s>\ 2 d„m 2 - 


2v 1 


v, v- 


|$| 6 

(3.15) 


leading again to ch = y 2 , and now to Ac6 = 3 y 2 9 m /v s . The parameter 9 m /v s = 5 
measures the relative importance of explicit vs spontaneous symmetry breaking in (3.1). 


Let us also comment on the case where the held s does not develop a vev and still 
the terms linear and cubic in s are absent from (3.1) (the symmetry G from (3.1) remains 
unbroken). In this case there is no Higgs-singlet mixing, but nevertheless the operator 
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|<L | 2 + |<I >| 2 is generated at 1-loop [57] with a Wilson coefficient ch = gg^ 2+2 (being 
n s the number of singlet scalar degrees of freedom). The universal suppresion of Higgs 
couplings can then be defined in a similar fashion to the previous case, y 1 = ch- 



+ obs/^SM 

ATLAS 77 7+8 TeV: /z gg F 

1.32 + 0.38 

ATLAS 77 7+8 TeV: +vbf 

0.8 ±0.7 

ATLAS WW* 7+8 TeV: /z ggF 

0.82 + 0.36 

ATLAS WW* 7+8 TeV: /x V bf 

1.66 + 0.79 

ATLAS ZZ* 7+8 TeV (inclusive) 

1 44+ 0 - 40 
i - 44 -0.33 

ATLAS bb 7+8 TeV: + V H 

0 2 +0 ' 7 
u ’ z -0.6 

ATLAS tt 7+8 TeV: /z gg F 

1 9+O.8 
1 * z -0.6 

ATLAS tt 7+8 TeV: +vbf 

1 6 +0 ' 6 
i - D -0.5 

CMS 77 (Mass Fit) 7+8 TeV: /z gg F+ttH 

1 iq+0.37 

1-1 —0.31 

CMS 77 (Mass Fit) 7+8 TeV: /xvbf+vh 

i ' 1D -0.58 

CMS WW* 7+8 TeV (0/1-jet) 

0 7A+V-' 2 ' z 
u - '^-0.20 

CMS WW* 7+8 TeV (2-jets, VBF tag) 

0 60 +0,57 
U.OU_ 0 .46 

CMS ZZ* 7+8 TeV (inclusive) 

0 oq+0.29 

U.yo_0 24 

CMS bb 7+8 TeV: /z V H 

1.0 ±0.5 

CMS rr 7+8 TeV (0-jet) 

0.34+1.09 

CMS tt 7+8 TeV (1-jet) 

1.07 + 0.46 

CMS tt 7+8 TeV (2-jets, VBF tag) 

0.94 + 0.41 


Table 2. ATLAS and CMS measured Higgs Signal Strengths /q = /z* bs //Zg M in h —» 77 [58, 59], 
h —> ZZ* [60, 61], h —> WW* [62, 63], h —> bb [64, 65] and /i 4 rr [ 66 , 67] final states. 


We now discuss the current experimental constraints on the mass of the singlet-like 
state and the parameter y, coming from the latest measurements of Higgs boson signal 
strengths from ATLAS and CMS, from oblique corrections to EWPO and from direct 
searches of heavy Higgs scalars at LHC. For the analysis of Higgs signal strengths, we 
consider the up-to-date measurements by ATLAS and CMS in h —> 77 [58, 59], h —> ZZ* 
[60, 61], h —> WW* [62, 63], h —> bb [64, 65] and /i -> rr [ 66 , 67] final states, shown in 
Table 2. We then perform a combined y 2 fit 


v2 _y 

Xh M A/Zj ) 


(3.16) 


with the various m = /^obs/^SM an d A /7 taken from Table 2, and with potential correlations 
among the different signal strength measurements not included in the fit. This yields an 
ATLAS and CMS combined limit y < 0.468 at 95%C.L. via Ay^(y) = Xhiv) ~ Amin = 
4, as shown in Figure 2 (horizontal solid-black). In addition, we obtain the 95%C.L. 
exclusion prospects for LHC at 14 TeV with C = 300 ftC 1 (horizontal dotted-black) and 
L = 3000 fb ” 1 (HL-LHC, horizontal dashed-black). In doing so, we assume that future 
measurements of Higgs signal strengths will yield //,- = 1 , and use the projected CMS 
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sensitivities 4 [68] AnXj n = 0.06, A//S 4 = 0.06, Artf™ = 0.07, A nZA = 0.08, Au™ n = 
0.11, A^ 0 = 0.40, A//3 qo 0 = 0.04, A/i^ = 0.04, A/ A || 0 = 0.04, A^ 00 = 0.05, 
A/4ooo = 0.05, A/i3Q 00 = 0.20. 



Figure 2. Present 95%C.L. exclusion limits in the (rn s , y) plane arising from ATLAS and CMS 
measurements of Higgs signal strengths as shown in Table 2 (horizontal solid-black) and from 
EWPO (blue). The shaded grey region is excluded at 95%C.L. by the combination (EWPO and 
Higgs signal strengths). The shaded yellow region may also be excluded by constraints from heavy 
scalar searches at LHC (dashed-dotted black), but these limits can be evaded in the presence of 
non-SM decays (see text). Also shown is the projected exclusion reach from Higgs signal strengths 
at the 14 TeV run of LHC with C = 300 fb _1 (horizontal dotted-black) and at HL-LHC with 
C = 3000 fb _1 (horizontal dashed-black), from measurements of the S and T oblique parameters 
with ILC-GigaZ and TLEP (blue) and their combination with the HL-LHC exclusion reach (red), 
and from a precision measurement of the Higgs production cross section in association with a Z 
boson a(hZ) (horizontal solid-green lines) at ILC (y/s = 250 GeV, C = 1150 fb -1 ) and TLEP 
(y/s = 240 GeV, £ = 500 fb" 1 ). 


Turning to EWPO, we perform a fit to the oblique parameters S, T, U using the best-fit 
values and standard deviations from the global analysis of the GFitter Group [69], with a 
SM reference point with rrq = 173 GeV and a 126 GeV Higgs mass. Under the assumption 
U = 0, this yields 

AS = S - S SM = 0.06 ± 0.09 _ / 1 0.91 

AT = T-T sm = 0.10 ±0.07 ~ y0.91 1 

4 These assume that theoretical uncertainties improve by a factor 1/2 compared to their present values, 
while all other systematic uncertainties improve by a factor 1 /y/C. 
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being ptj the covariance matrix in the S — T plane. The BSM corrections to S and T in 
the case of singlet-doublet mixing are given by 


AT = 


A S=-y z 

7T 

9 2 ..2 


ml 


16 7i Cw «em 


—Ht 


ml 


ml 


TTlr 




mt 


mi 


+ Ht 


ml 


m 


z/ J 


(3.18) 


with the functions Hg(x ) and Ht{x) defined in Appendix C of [70]. We then define 

A XEw( m s,y) = ^2 (A Oi(m s ,y) - A Of) (a 2 )^. 1 (A Oj(m s ,y) - AOj) , (3.19) 




where A O® denote the central values in (3.17) and (cr 2 )^- = GiPijGj , being Oi the S and 
T standard deviation from (3.17). We show in Figure 2 the 95%C.L. exclusion limit 
y(m s ) from Ax^ w (m s ,y) (blue), together with the 95%C.L. exclusion limit y(m s ) from 
the combination of A x\(y) and Ay| , w (m s ,y) (red). We also study the future exclusion 
reach that can be derived from prospects of measurements of EW precision observables 
by planned e + e~ colliders (see e.g. [71]): Assuming a SM best-fit value, the ILC GigaZ 
program’s expected precision is as = 0.017 and gt = 0.022 [69], while measurement of 
EWPO at TLEP could yield as = 0.007 and gt = 0.004 [72, 73]. Figure 2 includes the 
95% C.L. exclusion reach for y{m s ) both for ILC and TLEP (blue), as well as the respective 
95%C.L. exclusion reach when combined with HL-LHC (C = 3000 fb” 1 ). Figure 2 also 
shows the potential 95% C.L. exclusion reach from a precise measurement of the Higgs 
production cross section in association with a Z boson a(hZ) at ILC with y / s = 250 GeV 
and C = 1150 fb _1 (with Aa(hZ)/a{hZ) ~ 0.012) [74] and at TLEP with -^/s = 240 
GeV and C = 500 fb _1 (with Aa(hZ)/a(hZ) ~ 0.004) [72], From Figure 2 it is evident 
that this precision measurement is potentially the most powerful probe of a non-SM Higgs 
admixture, significantly surpassing the reach of EWPO measurements. 

It is worth stressing that the contributions to the S and T oblique parameters in the 
Higgs portal scenario could be re-derived in the EFT approach via operator mixing under 
renormalization group evolution (RGE), as shown e.g. in [75]. Specifically, the running of 
ch from the matching scale m s down to mz generates a contribution to S and T, given at 
leading order by [75, 76] 

AS = -*-Ctf(m s )log (—\ AT = - 3 g c H {m s ) log (—) (3.20) 

6 7 t {m z ) 87 rc^ \m z ) 

which correspond to the leading order contributions in (3.18) for m s 3> mz- 

Finally, Figure 2 includes the latest 95% C.L. limits on y(m s ) (dotted-dashed black 
line and yellow region) obtained from CMS searches of a heavy neutral scalar decaying to 
ZZ —> 4£, ZZ —» 2£2j and ZZ —> 2£2u final states [77] (searches in decays to WW —> 2£2u 
final states [78] are found not to be as sensitive). These constitute at present the most 
stringent constraint on y for m s < 500 GeV, but we stress that for m s > 2 mh = 250 GeV 
these limits may be weakened/avoided for a significant branching fraction Br(/i 2 —> h\h\). 
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Figure 3. Present and future (as defined in Figure 2) 95% C.L. exclusion limits on the size of 
ch{itiz) (Left) and Xce(mz) (Right) as a function of m a , from the combination of Higgs signal 
strengths and EW oblique parameters. For the case of Xce(mz), the region leading to a strongly 
first order EW phase transition as required by EW baryogenesis is shown as a horizontal orange 
band, while the grey (brown) region is currently excluded at 95% C.L. for S = /r m /u s = 1 (6 = 5). 


From the present limits and future prospects discussed above, we can derive cur¬ 
rent/projected bounds on the size of the Wilson coefficients c/y(mz) and c§{m z ) in this 
scenario. From ch(iti s ) = y 2 and Ac6(m s ) = 3 y 2 fim/v s , we perform a 1-loop RGE [28] (see 
also [29]) 


CH(m s ) 


CH(m z ) - 



2 


3 g ' 2 + 24A + 12 y 2 \ 

16 7r 2 


c H (rn s ) log 



(3.21) 


A c 6 (m s ) 


[_ 27 0 2 

A c e (m z ) - -—-— 


\g ’ 2 + 108A + 18 y 2 ] 

16 7r 2 

[-3 g 2 + 40 A] 
+ 8vr 2 


A c 6 (m s ) log 
Xc H {rn s ) log 



(3.22) 

(3.23) 


Current bounds on CH(m z ) and ce(m^) as a function of the new physics scale m s , as well as 
the projected 95% C.L. exclusion sensitivity from the combination of HL-LHC and GigaZ, 
and HL-LHC and TLEP EWPO ( TeraZ ) measurements, are shown in Figure 3. 

As recently noted in [75], a sizeable CQ(m z ) could in this scenario lead to a strongly 
first order EW phase transition, potentially allowing for EW baryogenesis in the Early 
Universe [79-81]. Requiring a sufficiently strong first order phase transition imposes a 
lower bound Ace(m^) > 0.105 5 , while the successful completion of the phase transition 

5 The EW phase transition occurs at T = T c ~ 100 GeV, and we can approximate T c ~ mz for cq. 
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sets the upper bound Xce(mz) < 0.211 [80]. As shown in Figure 3 (Right), the present 
combination of measured Higgs signal strengths and EW precision data rules out a strong 
EW phase transition, except for very large values 6 = y m /v s 3> 1 combined with relatively 
low singlet masses m s . These results are somewhat stronger than those of [75], the reason 
being that the analysis of the EW phase transition in the presence of cq from [79] (with 
reference [75] follows) considers thermal masses as the only thermal effects in the effective 
potential, and tends to overestimate the strength of the EW phase transition as compared 
to a full 1-loop analysis [80]. Finally, while we stress that our results certainly do not rule 
out a strong EW phase transition in the singlet scalar extension of the SM (see e.g. [45]), 
they imply that having it originate from the presence of cq is currently challenging, and 
the entire region of parameter space will be covered by the combination of HL-LHC and 
GigaZ. 

Let us finish this section on the extension of the SM Higgs sector via a singlet scalar by 
summarizing the results. We have seen that the EFT leads to just two non-zero operators 
at tree-level, 


c H = y 2 (mixing) 

n s (X m v\ 2 ■ • \ 

CH = wp- UrJ (no mixmg) 

A ce = 3 Sch (only w. explicit symmetry breaking) 


(3.24) 


3.2 Two Higgs Doublet Models 

Theories with two scalar doublets appear in a wide variety of scenarios, ranging from the 
MSSM and its extensions [82, 83] to models of dark matter [84-86] and neutrino masses 
(see e.g. [87]). The phenomenology of Two Higgs Doublet Model (2HDM) scenarios has 
been widely studied in the literature (for a general review of 2HDMs, see [ 88 ]), including 
the impact of measured properties of the observed Higgs boson on the 2HDM parameter 
space [89-94, 96, 97], the possible connection to the EW phase transition and baryogenesis 
[98-101] and its potential LHC signatures [102-108]. 

We consider here 2HDM scenarios with a ^-symmetry which is at most softly broken, 
avoiding tree-level Flavor Changing Neutral Currents (FCNCs) in the Yukawa sector [109]. 
This ^ 2 -symmetry leads to four types of 2HDM (see e.g. [ 88 ] for details), according to the 
way in which T 12 are coupled to the different SM fermion species, as shown in Table 3. 
The scalar potential for the two scalar doublets <Li ! 2 reads 


Ai| ,4 
+ — $i 
2 1 1 


H tree ($i, $ 2 ) = n\ |$i| 2 + l4 l$2| 2 - h 2 $I$2 + h-c. 

A 5 (V^y + h.c. 


+ ^|4 > 2| 4 + A 3 |4>i| 2 |4> 2 | 2 + A 4 <L t 1 4>2 


1 

+ 2 


(3.25) 


In the following, we consider CP-conserving scenarios and set A 5 and y? to be real. The 
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Up-type Quarks 

Down-type Quarks 

Charged Leptons 

Type - I 

$2 

( I >2 

$2 

Type - II 

$2 

$1 

$1 

i — Specific 

4>2 

<f >2 

4»i 

Flipped 

4*2 

$1 

4?2 


Table 3. Classification of 2HDM with a Z^-symmetry in the Yukawa sector. For each 2HDM-type, 
we indicate which scalar doublet couples to up-type quarks, down-type quarks and charged leptons. 


doublets and their vevs at the EW minimum can be written as 

= ( h fX lk ) , ($i) = (° ) , <$ 2 ) = (°), (3.26) 

V C2 / VV2/ \V2j 

with v\ + v\ = v = 246 GeV and v 2 /vi = tan/3. The angle /3 also parametrizes the rotation 
to the mass eigenbasis for the charged states and neutral CP-odd G°,A° states, 

with G ± ,G° being the Goldstone bosons and H±,A® the physical states. We likewise 
define a to be the mixing angle parametrizing the rotation to the mass eigenbasis for the 
CP-even neutral states h,H° (see Appendix A.l). The h parameter is responsible for the 
soft-breaking of the ^-symmetry in (3.25). After EW symmetry breaking, the parameters 
Hi, A i in (3.25) may be written in terms of the masses of the physical states rrih , m H o, m^o, 
rriff ±, the mixing angles a , (3 and the /i parameter, as shown in Appendix A.l. 


In order to obtain an EFT for the SM Higgs, we may perform an SU{ 2) rotation 
from the field basis < hi, < J >2 to a basis H\. H 2 in which (Hi) = and ( H 2 ) = 0 (the so- 
called Higgs basis). This rotation is precisely parametrized by the angle /3. After the field 
rotation, the scalar potential for H\ 2 reads 


V lree (Hi, H 2 ) = h\ \Hi\ 2 + ~h\ \H 2 \ 2 - h 2 H\H 2 + H.c. 


, I TT | 4 


A 2 1 TT ,4 , S: I tt |2 | tt |2 , t trtrr 2 , ^5 


+ y \H 2 \* + A 3 |fTi|' | H 2 \ z + A 4 H\H 2 


+A6 


\H\ | 2 H\H 2 + H.c. 


+ T 

+ A7 


h\h 2 ) + H.c. 


\H 2 \ 2 h\h 2 + H.c. 


(3.27) 


with fli, fl and Aj being functions of the original parameters in (3.25). We note that the 
field rotation may generate A 6,7 even if initially absent from (3.25) due to the Z 2 -symmetry. 


We may then construct an effective Lagrangian for H\ by matching to the theory with 
the second doublet H 2 at the scale /I 2 , assuming jl 2 i$> v 2 . We stress however that the 
doublet Hi can only be fully identified with the SM Higgs doublet in the alignment limit 
cos (/3 — a) = cg_ Q —> 0, where mixing in the CP-even sector is absent (see the discussion in 
Appendix A.l). Away from alignment (a 7 ^ j3 — 7r/2), these mixing effects lead to tree-level 
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modifications of SM Higgs couplings (see e.g. the discussion in [110]). For a general 2HDM, 
C/ 3 - a = 0 is possible at tree-level 6 [111, 112], and indeed the latest analyses of Higgs data 
by ATLAS and CMS strongly prefer c 2 _ Q <C 1, as shown in section 3.2.3. 

In the following, we introduce and discuss two benchmark scenarios for the use of 
Higgs Effective Theory in 2HDM: An exact alignment scenario C/ 3 - a = 0 and a scenario 
with |cg_ Q | -C 1 (MSSM-like scenario). We then analyze the constraints on the 2HDM 
parameter space obtained from measurement of Higgs signal strengths and EW Oblique 
Parameters, and the corresponding constraint on the values of the Wilson coefficients of 
the Higgs EFT. 

3.2.1 Benchmark A: Exact Alignment C/ 3 - a = 0 

As discussed above, in this scenario mixing effects are absent and Hi is precisely the 
SM Higgs doublet <f>. The scalar potential (3.27) simplifies in this limit, since A@ and p 
are both oc cg_ Q and vanish as shown in (A.16)-(A.17). We then match the 2HDM to an 
SU{2)l x U(l)y invariant EFT with the field content of the SM, being H 2 the only massive 
field which will be decoupled in the matching calculation. The D = 6 effective operators 
from (2.1) are generated first at 1-loop order, with the corresponding Wilson coefficients 
given by 


ch = ~ 


— 4 A 3 A 4 + A 4 + Ag — 4 A 3 


. 192 7 r 2 p\ 


cq (A4 T A 


192 7T 2 p 2 

„,2 


CT = (A 1 - Ag 


192 7 r 2 p 2 


mw A 3 
2567r 2 p 2 


rn 


cw = —chw = 


w (2 A 3 + A 4 ) _ 8 _ 

A ^ 2 ~ 3 C7 + 1 Q9 -tt 2 r,2 


'w'w A 4 


cb = -CHB = 
C2W 


_ m w 


192 7T 2 /i 2 
(—2 A 3 + A 4 ) 


192 7 r 2 


192 7r 2 p4 


= _8 rn 2 w A 4 
3 7 192 7T 2 P 2 


m 


C3W = 


w 


14407r 2 p z 2 


from which we immediately obtain 


cw + cb = 


’cn'w 


16 


mL A 3 


c w - cb = — Cry = 2 

3 48 7T 2 H 2 


(3.28) 


(3.29) 


96 7T 2 p 2 

Before we continue, let us note that Cy, cw and cb may take positive and negative values, 
as the bounded-from-below conditions (A.2) do not restrict either possibility. The various 


b In contrast, for the MSSM c/ 3 - 


0 is only obtained at tree-level in the decoupling limit P 2 °°- 


Moreover, cp- a = 0 at loop level is only possible in a very small portion of parameter space [111, 112], 


- 16 - 


















relations among the Wilson coefficients chw, cw, chb, cb and c 7 in (3.28) imply that the 
D = 6 effective operators Ohw > Ow, Ohb, Ob and 0 7 can in fact be re-casted in terms 
of three operators 

W*,WT + [H\T 2k H,]W^ + E, 0, t (3.30) 

in w L w 

with 

^-|Fi| 2 wl~;jvl w = O ww = 4 (Ow - Ob + O hb - O hw ) + 0 1 (3.31) 

171 w 

2j] 4- [h\T 2k H-i] = O wb = 4 (Ob - O hb ) - 0 7 (3.32) 

m w 

The three operators in (3.30) share the common feature that they do not involve derivatives 
of the Higgs field Hi. This very property is in fact responsible for the relations among 
Wilson coefficients in (3.28), since the Feynman diagrams involved in the EFT matching 
(see Figure 4-Left) do not involve Higgs field derivatives. 


3.2.2 Benchmark B: Departure from Alignment cp- a <C 1 (MSSM-Like) 

Upon departure from the alignment limit, mixing in the CP-even sector leads to several 
effects that are absent for cp_ a = 0: First, there is a modification of the couplings of 
the Higgs boson h to gauge bosons and fermions at tree-level, parametrized in terms of 
k -factors for vector bosons Ky = 9 kv^v v /and fermions Kf = 9hff/9hjf (f° r u P~type 
quarks k u , down-type quarks Kd and leptons Kg). For the various Types of 2HDM from 
Table 3, these are given in terms of cp- a and tp = tan/3 as 


Type — I : Ky 
Type — II : Ky 
£ — Specific : Ky 
Flipped : Ky 


_ _ _ C |S—a . 

S/3—a j fi'U ^d k i , V S/3—a 

t/3 

_ c fi—ot . _ _ 

S/3—a j ^u , T SR— a , Kd Kg SR—a 

t/3 

_ . _ c /3 ~‘a | . _ 

S/3—a j — , T SR—a j — SR— a 

tp 

_ _ c /3 —a . _ 

Sj3—a i £3, u — Kg — — h Sp—a , Kd — Sp— a 

tp 


tp Cp—ot 
tp Cp—a 
tp Cp—a 


(3.33) 

(3.34) 

(3.35) 

(3.36) 


In addition, away from alignment g B ohh (coupling between H° and two Higgs bosons h), 
9 h°v V v (coupling between H° and two gauge bosons V),) and gpy^h (couplings among h, a 
gauge boson V p and (j) = A 0 , H±) are non-zero, giving rise to 1-loop diagrams contributing 
to the vertices V p V u h and V P V U V P with both heavy (H°, A 0 , iL 1 * 1 ) and light (h, V p ) states 
running in the loop. 


Let us discuss all these effects in an SU(2) B x U( l)y invariant Higgs EFT approach. 
We first note that (H\) = and (H 2 ) = 0 in (3.27) imply the relations 


Mi = 


Xiv 2 

2 


~2 \ 2 
M = A 6 v 


(3.37) 
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which may be simply regarded as minimization conditions. The mass matrix for the neutral 
CP-even states is non-diagonal, brought into a diagonal form via the rotation matrix U 


U = 



p 2 3 Ag v 2 Ag v 2 
2 A2 2 Pl p2 


(3.38) 


The deviation from alignment is then parametrized by (Xqv 2 ) / p2, which matches the cor¬ 
responding expression for cp_ a after EW symmetry breaking (see Appendix A.l) 


r 2 

C/3—a 


2\2 


{XqVj 

~A 


V4 < 1, 

M 2 


(3.39) 


recovering the well-known scaling result for c|_ a in 2HDMs as the decoupling limit is 
approached [111]. 

At tree-level, the EFT matching generates effective operators suppressed at least by 
1 / p 2 (except for eg, which receives a further contribution suppressed only by I/P 2 ), with 
operators of different D contributing to order v^/p^- 

£eF 3 ^-{D h\){ D ^ H i ) + c ^-d^\H 1 \ 2 d fl \H 1 \ 2 
M 2 M 2 

+ ^D^HIIHP^D^HPHP 2 ) - (hf ~ |#i| 6 • (3.40) 

M 2 \ M 2 ^ M 2 / 


The matching procedure yields cg^ = A§, c^h = ch = 0 and C 2 h = 1, and use of the relation 
p 2 = Ag v 2 results in both the first and the third term in (3.40) contributing at order 
(Xqv 2 ) 2 / p2- The first term rescales both the SM Higgs kinetic term and its couplings to 
gauge boson by the same amount, so it does not have a net effect. The third term does 
however include a rescaling of the SM Higgs kinetic term oc c|_ a that is not compensated by 
a similar one in the gauge boson interactions, leading to the well-known tree-level deviation 
from the SM Higgs couplings to gauge bosons away from alignment, proportional to c^_ Q . 

Regarding eg, the presence of a non-zero Ag also yields an extra contribution at 1-loop 
as compared to (3.28). The full result for eg away from alignment is 


(A| + A| + 12 Ag) v 2 
192 7 r 2 P 2 


X 2 v 2 Ag p 2 v 2 \ 

Am! m! J 


(3.41) 


Let us also stress that similarly to the tree-level effects just discussed, the extra 1- 
loop diagrams appearing away from alignment (involving both light and heavy degrees of 
freedom) are proportional to c 2 p_ a - Thus the contribution of these diagrams to the various 
Wilson coefficients from (2.2) and (2.5) is at least of order (v^/p^) x 1-loop and can be 
safely neglected. As a result, the values of the Wilson coefficients ct, chw-, cw, chb , cb, 
Cry for c|_ ct <C 1 remains unchanged w.r.t. the alignment scenario. 

The interactions between the SM Higgs and fermions receive tree-level corrections of 
the form shown in Figure 4-Right away from the alignment limit, encoded in the D = 6 
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Figure 4. Left: Feynman diagram responsible for anomalous Higgs couplings to gauge boson at 1- 
loop. Right: Feynman diagram responsible for anomalous Higgs couplings to fermions at tree-level. 


effective operators 


C 


( 6 ) _ 

SILH — 


with 


C ^H\H l Q l h\u r + ^ H\H r Q L H l d R + C ^h\h ] L l H x l R 

L V z V z V z 


- n . _ m u fu X e v 2 _ __ _ m d f d X e v 2 _ _ _ mefi X 6 v 2 

Cu Vu ~ o C-d Ud ~ 9 (: f IH 




V v « 

and f u d i depending of the 2HDM-Type under consideration, and given by 


Type - I : f u = f d = f e = tp 1 
Type -II : f u = t ^ ; ,f d = ft = ~tp 
t - Specific : f u = f d = t^ 1 ; f e = -tp 
Flipped : f u = fe = ; f d = -tp 


These results reproduce the /{-factors from (3.33-3.36) to leading order in cp- a . 


(3.42) 


(3.43) 


(3.44) 

(3.45) 

(3.46) 

(3.47) 


3.2.3 2HDM Constraints: Higgs Signal Strengths & EW Oblique Parameters 

The ATLAS and CMS measurements of Higgs boson signal strengths from Table 2 constrain 
the allowed parameter region of the 2HDM in the plane (cp_ a , tan (3) (see e.g. the analyses 
from [89, 91-93, 96, 103]) through a combined x 2 fit, which in this case treats the different 
signal strength measurements as independent (not correlated) 

= (3.48) 

Each jli (the expected signal strength for each channel) may be expressed in terms of 
rescaling k- factors for the coupling of the Higgs boson to vector bosons Ky, up-type quarks 
k u , down-type quarks K d and leptons K£. The expression of each Jli corresponding to each 
process considered in Table 2 in terms of these /{-factors is given in table 4 
with Kj defined in (3.33-3.36) and k 7 (k u , Ky) = 1.26 Ky — 0.26k u , and T(Kd, Ky, k u , k#) 
given by 

T(/e d , k v, n u , Ki) = 0.2427 Ky + 0.1124 /c 2 + 0.578 K 2 d + 0.0637 k\ (3.49) 


19 


















P = p(k)/psm 

IT TggF 

X K 7 (fv u ,fvy ) 2 X T (n d , KVTuTl)~ l 

77 : /avBF 

Ky x Kj(k u ,k v ) 2 x T(K d: nvTuTi)~ l 

WW*: /UggF-O/l-jet 

^“u X Hy X r(K d ,Ky,K u ,K£) _1 

WW*: p v bf 

Hy X r(K d ,Ky,K u ,K£) _1 

ZZ*: Inclusive 

kI X X r^Ky,^,^) -1 

bb: /i VH 

Ky X K 2 d X r (K d , KV, K u , Hl)~ X 

rr: /i gg F-0/l-jet 

H 2 u X /cf X V(K d , K V , K u , Kt )- 1 

tt : pvbf 

Hy X K 2 X r(K d , Ky, K n , Kf) 1 


Table 4. Modified (by K-factors) Higgs Signal Strengths p t . 


We present the 68 %, 95% and 99%C.L. limits on the (cp- a ,tp) plane from global fits 
of the light Higgs boson couplings, given respectively by A x 2 h (cp- a ,tp) = x 2 h (cp~ a ,tp) - 
Xmin = 2 - 27 > A Xh( c P-a,tp) = 5.99 and A xl( c p-a,tp) = 9.21, for the different Types of 
2HDM in Figure 5. These are particularly stringent for Type II and Flipped as compared to 
Lepton-Specific and specially to Type I, due to the different nature of the Higgs coupling to 
down-type quarks. Figure 5 also includes the projected 95% C.L. exclusion limits for LHC 
at 14 TeV with C = 300 fb -1 (dashed-black) and C = 3000 fb _1 (HL-LHC, dotted-black) 
assuming the measured Higgs signal strengths being m = 1 and using the projected CMS 
sensitivities from Section 3.1. 

We now analyze the constraints from EWPO, performing a fit to the oblique parameters 
S, T under the assumption U = 0 using the best-fit values and standard deviations from 
the global analysis of the GFitter Group [69] (m* = 173 GeV and = 126 GeV), as 
shown in (3.17). The 2HDM contributions to A S and AT are given by [113] 


AS 


a 2 s 2 
y b w 


96 7r 2 «em 


+ c 


(3—a 
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(3.51) 


with G a> b,c, G A) b and F A b given in Appendix A. 2 . We note that both AS and AT are 
independent of tp. We define 


m 2 o m 2 0 

Xo = o , x A = —f- , 
m H± m H± 

and construct a A'x 2 EW function as 


(3.52) 


Ax 2 EW (m 2 H ±,xo,x A ,cp- a ) = J2 (A Oi - A 0°) (<r 2 )r. 1 (A Oj - AOj) , 


(3.53) 
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Figure 5. 68 % C.L. (yellow), 95% C.L. (green) and 99% C.L. (red) allowed regions in the (cos(/3 — 
a), tan/3) plane from up-to-date measurements of Higgs boson couplings (see text for details), shown 
respectively for 2HDM of Type I (Upper-Left), Type II (Upper-Right), Lepton-Specific (Lower-Left) 
and Flipped (Lower-Right). In each case, the projected 95% C.L. exclusion limit for LHC at 14 
TeV with C = 300 fb -1 (dashed-black) and C = 3000 fb” 1 (dotted-black), assuming the measured 
Higgs signal strengths being ji t = 1, is also shown. 


with A Of and (u 2 )jj defined after (3.19). The results of the fit are shown in Figure 6 for 
various values of cg_ a and m H ± = 400,600 GeV. 

The calculation of the oblique S and T parameters in the 2HDM matches with the 
D = 6 effective operators present in the EFT. Considering first the alignment limit and 
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Figure 6. 95% C.L. allowed region from the 2HDM fit to S,T oblique parameters (3.53) in the 
plane (xq, xa) for mn± = 400 GeV (Left) and mu± = 600 GeV (Right). The allowed regions 
correspond respectively to cp_ a = 0 (grey, light), cg_„ = 0.2 (red, dark) and cg_ a = 0.4 (orange, 
medium). 


expanding (3.50) and (3.51) around xq,xa ~ 1 (as indicated by the present 95% C.L. 
experimental limits), AS and AT in the 2HDM read 


AS 
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(3.54) 


where we note that the functions Ga,b,c scale as (1 — xa,o) 2 and can be neglected with 
respect to the log(xo /x A ) term in AS. By means of (A.14) and (A.15) in Appendix A.l, 
the above relations can be rewritten as 
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(3.55) 


For m 2 H ± u 2 , and upon the replacement rn 2 H ± —> fl 2 , these match the EFT results (3.28) 
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where we have also translated into the Altarelli-Barbieri parametrization [114] for reference. 
The bounds on these effective operators are [69] 

c T {mz) € [~1.5,2.2]xl0- 3 and (c-wipiz) + c B (mz)^J G [—1-4,1.9] x 1CT 3 . (3.57) 

Away from alignment, AS and AT receive further contributions oc c^_ a . In an EFT 
language, these are captured via RG evolution: various operators from (3.40) may mix 
with Ow, O b and Ot via RG running, as is the case of e.g. Oh [28], whose mixing at 
leading order would yield 

c T (m z ) - cr{fa) - c H {fa) log (—) (3.58) 

32 -k z \mz J 

cw(rnz) + c B (m z ) - c w (jj, 2 ) + c B {fi 2 ) + c H {fi 2 ) log ( — ) . (3.59) 

24 7r z \mz J 

However, for the 2HDM ch{^- 2 ) = 0 at tree-level, and thus the Wilson coefficient responsible 
for the contribution to AS and AT dependant on c|_ a is c^h [27]. 


Electroweak+LHC constraints 



Figure 7. 95% C.L. limits on the 2HDM parameters A 4 and p. 2 - The red-dotted lines enclose the 
region allowed by constraints on the S parameter, with the black-dashed lines the corresponding 
LHC limits from [21]. The purple-solid region corresponds to the combined limits. 


3.2.4 Alignment limit: complementarity of EWPTs and LHC limits 

In the previous section we discussed how LEP electroweak constraints translate into the 
full theory (3.28), and in turn in the effective theory (3.56) which affected the operators 
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c-t and a combination c\y + c-b- It is time to move onto the constraints from LHC, the 
most stringent ones on operators affecting the decay of the Higgs to photons (c 7 ) and 
production through gluon fusion ( c g ). Other operators are better determined by looking 
at other production and decay channels. In particular, three combinations [19] 

cw — cb , chw and chb , (3.60) 

are constrained with the help of Higgs production in association with a vector boson, and 
with diboson data [21], In the 2HDM, the global fit is more constraining than a fit for a 
general EFT at dimension-six level. In particular, by varying all the free parameters in the 
2HDM simultaneously in a fit to Higgs plus diboson data, one obtains the following 95% 
C.L. regions [21] 


chw £ (0.0004,0.02) 

c g <E -(0.00004,0.000003) 

c 7 G (-0.0006,0.00003) (3.61) 

Whereas the constraint on the T parameter favours a A 5 close to A 4 , the constraint 

on c 7 from the LHC restricts the values of A 3 . The most interesting case is when we 
compare constraints from LEP and LHC which set limits on the coupling A 4 and the scale 
of new physics. Indeed, one can see this interplay between the two sources of data in 
Figure 7, where we present results in the A 4 and /I 2 plane. The red-dotted lines enclose 
the region allowed by constraints on the S parameter, and the black-dashed lines are the 
corresponding LHC limits from [21]. The purple-solid region corresponds to the combined 
limits and it is quite constraining as the two regions have little overlap. Note that the 
preferred region at 95% C.L. is consistent with the presence of light particles. This is due 
to the slightly non-zero positive value of chw hr the global LHC fit and the tension with 
LEP limits, which favour negative values of chw ■ This hint for new physics is gone at the 
level of 3cr. 

Finally, we comment on the 2HDM and the impact on the EW phase transition. As 
recently shown in [106], a first order EW phase transition strongly favours a large 
splitting, which via (A. 17) implies A 5 < 0 and sizable. Since this implies that it is not 
possible to simultaneously have vanishing AS and AT unless jj ,2 —> 00 , it would be possible 
to observe a deviation from the SM in AS or AT as the experimental precision increases 
(see Figure 8 ). 

3.3 A Dilaton/Radion Scenario 

So far we have discussed scenarios where the operators involving the Higgs and massive 
vector bosons, such as c\y, are generated at loop level. In this section we present a case 
where they appear at tree-level, through the exchange of a dilaton or a radion scalar 
particle. We consider two equivalent scenarios: 

• The extra-dimensional radion: A radion r is the excitation of the graviton 
perturbation in extra-dimensions along the direction of the extra-dimension. The 
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interactions of the radion with SM particles are then obtained by expanding the 
metric at linear order in r. We will consider conformal metrics of the form 


ds 2 = w(z) 2 (r]^ l/ dx fJ, dx l/ — dz 2 ) 


(3.62) 


where w(z) is the warp factor, which depends on the extra-dimension z. Minkowski 
space corresponds to w = 1 and in Anti-deSitter metrics w = 1 /(kz), with k the cur¬ 
vature in the extra-dimension. The extra-dimension is compactified, with z stretching 
between two points which we will call the UV and IR branes, i.e. z G [zuv■, zir\ . The 
interaction then reads 

S = f d d x\/—gC D f d d Xy/^gw 2 (z)2rT^ (3.63) 


where T'ft is the trace of the stress-energy tensor. The dimensional reduction of the in¬ 
teraction in (3.63) to a four-dimensional (4D) interaction depends on the localization 
of the radion and the SM particles, and generally speaking it would look as 

C iD =^rT (3.64, 

where now r and T l = Tr(T^ u ) are the dimensionally-reduced radion and trace of the 
stress-energy tensor of species i. A corresponds to the compactification scale, and 
the coefficients q account for the overlap of different bulk species i with the radion. 

In warped extra-dimensions, the radion is a field localized near the IR brane at zjr. 
Therefore, for fields localized near the IR brane, c* ~ 0(1), for fields delocalized in 
the bulk (with a flat profile) one finds c* ~ rz uv and for fields on the UV 

brane, c* with a a positive number [115]. As an example, in Anti-deSitter 

(AdS) models, the suppression for bulk fields is f w(z)dz = log yjy, which in Randall- 
Sundrum translates into log(Mp/TeV) ~ 0(30). 

Although massless gauge fields do not contribute to the trace of the stress-energy 
tensor at tree-level, loop-level couplings are generated by the trace anomaly 

Tv = — V s ^ lCXl U* f3 65i 

H,anom / j ) 

i 


where i runs over the gauge groups of the SM, SU( 3) c x SU(2)l x U(1)y, on = g 2 /4tt 
and gi are the respective coupling constants. The bi are the /3 function coefficients 
leading to anomalous scale-invariance violations (dilaton) or, equivalently, contribu¬ 
tions from fields localized near the IR brane or bulk running (radion) [116, 117]. The 
values of bi depend on the degree of compositeness of fermions in the SM and possible 
new CFT contributions. For simplicity we neglect these and use the pure CFT value, 
bf FT = 8 tt 2 /(gf log(zuv /zir)). In this limit the coefficient of the anomaly is then 
independent of the value of the gauge coupling 


rpil 

x H,anom 


log {zir/zuv) 




Z7* 


(3.66) 
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• The dilaton of spontaneous breaking of scale invariance: Assume the exis¬ 
tence of a new sector whose couplings are scale invariant, but this scale invariance 
is spontaneously broken by some dynamics. For example, in QCD scale invariance 
is spontaneously broken by confinement dynamics, leading to quark condensates. In 
our example, we are considering a new strongly coupled sector, possibly with similar 
dynamics as QCD. The Goldstone boson of this spontaneous breaking is a dilaton. 
Let us denote the dilaton by r as in the radion case. Non-derivative couplings of the 
dilaton to SM fields are proportional to the explicit breaking of conformal invariance, 
namely they are of the form 


jd^ (3.67) 

with / the scale of the spontaneous breaking of dilatation symmetry and J^ the 
global conserved current for this symmetry, J^ = T^ v x v . Therefore the coupling of 
the dilaton to the SM fields take the same form as those of the radion (3.64) [118]. 
We stress that the quadratic couplings to SM particles do not have the same form 
for the radion and dilaton [118], but this difference will not affect the matching of 
dimension-six coefficients. 


The radion/dilaton mass is a model-dependent parameter, related to the mechanism 
of stabilization of the extra-dimension or the explicit breaking of dilatation symmetry. In 
absence of stabilization the radion is massless but one could stabilize this dilaton/radion 
in several ways, and no definite prediction of the mass can be drawn unless we specify 
the mechanism. For example, in the Goldberger-Wise mechanism [119] in warped extra¬ 
dimensions the radion mass is a function of the vev and mass of the stabilizing field [120] 
and could be very light as well as around the scale of compactification (or spontaneous 
breaking) /. 

We proceed to integrate out the dilaton assuming its mass m r is larger than the scales 
we are probing with colliders. The effective Lagrangian has the form 

c ‘" = ~T ^ t2 <3 ' 68) 

where T is the trace of the stress-energy tensor. For the Higgs and gauge bosons, 

T C -2\Dn<S>\ 2 + 4H(4> t 4>) - — (3.69) 

87T m 

where V = —m 2 |<]?| 2 -F A|<I>| 4 . One can then extract the Wilson coefficients of the effective 
operators 
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(3.70) 


We stress again that for a pure CFT value bf FT , b\oti = 6202 = b^a^. We could also 
consider a more general situation, where the coupling of the dilaton/radion is not universal. 
This amounts to setting different values of the coefficients q in (3.64) and would lead to 
a prefactor in the coefficients of the effective operators (3.70), dependent of the degree 
of overlap of the wavefunctions in the bulk (radion) or participation on the composite 
dynamics of the species (dilaton). 


While the contributions to AT vanish in this case at tree-level [121, 122], they will be 
generated both at 1-loop (this is also the case for AS) and through operator mixing due 
to RG running. The 1-loop contribution to AT and AS is given by [121] 
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(3.71) 


The corresponding contributions to cximz) and cw(jnz) J rc.B{jnz) from the RG running of 
ch are given by (3.58-3.59) with the substitution p ,2 -A m r and using c#(m r ) from (3.70). 
We note that the RG running gives a contribution to AT and AS which despite begin 
suppressed by m 2 h /wi 2 w.r.t the 1-loop contribution (3.71), may become dominant due to 
a much smaller numerical suppression. The fact that the two contributions have opposite 
signs leads to a partial cancellation effect. In Figure 8 we show the results on the AS — AT 
plane, and compare the correlation in the oblique parameters with that appearing in the 
ones encountered respectively in the Higgs portal and 2HDM scenarios from Section 3.1 
and 3.2. 


4 The validity of the EFT: from operators to total rates and distributions 

The EFT approach, where higher-dimensional operators are written as a way to encode 
effects of New Physics, is a good approximation at low energies and should not be used 
in an arbitrary range of energies. Indeed, as these operators are suppressed by some New 
Physics scale, parton level cross-sections would diverge with the parton energy. 

This growth with energy is a signature that the EFT approach breaks down when 
one is able to probe the dynamics of the heavy particles one integrates out. In hadron 
colliders, this corresponds to the moment when the partonic energy y/1 is comparable with 
the masses in the loop. Specifically, for a tree-level exchange as in Figure 1, 

Vi — M (4.1) 
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Figure 8. AS — AT constraints on the radion/dilaton (black), singlet-Higgs portal (green) and 
2 HDM w. custodial symmetry (brown) scenarios. The thick (thin) ellipses correspond to the 
68% C.L. ( 95 % C.L.) allowed region for LEP (blue) and the projected ILC Giga-Z (yellow) and 
TLEP TeraZ (red), see Section 3.1 for more details. The zoom on the percent region (Bottom) 
includes points for the Higgs portal with various values of mixing parameter y and mass m s , for 
the radion/dilaton with / = 2 TeV and various values of the radion mass m r , and for the 2 HDM 
with A4 = A5 = —7r and several values of fi.2- 


The same argument applies to the dilaton exchange, and to fermion couplings in 2HDM 
as in Figure 4. On the other hand, in loop-induced processes as in Figure 10 the validity 
of the EFT extends to the threshold to produce a pair of new particles, namely 

~ 2 M . (4.2) 
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Figure 9. Higgs produced in association with a vector boson: Invariant mass distribution. Grey and 
pink distributions correspond to SM and EFT with chw = ~cw = 0.004. The width correspond to 
varying the factorization scale /r = rrih + mz in a range of 2/r and /i/2. The value of c\y is consistent 
with Runl data and will be explored in the Run2 LHC. The value 6 is the deviation respect to the 
SM central value. 


How this energy is related to the strength of the dimension-six deviation from the 
SM depends on the model. Let us focus on operators of the type cyy and chw which 
contain derivatives on the Higgs field. We have found the matching with the 2HDM and 
radion-exchange is roughly speaking 

~ ^ 9^2 Jp ( 2HDM ), and ^ ^ (radion/dilaton), (4.3) 

where M plays the role of the mass of the extra Higgses and dilaton, A denotes a combi¬ 
nation of quartic couplings in the Higgs potential, and A an effective scale suppression for 
the dimension-five interaction between the dilaton and the SM particles, A > M. 

At energies around the mass threshold , the EFT description should be substituted by 
the full UV theory including the resonance dynamics. The question we would like to ask 
is what is the deviation of the EFT predictions from the full UV model, which depends on 
the process and distribution one is looking at. A full comparison of the UV theory and the 
EFT results is beyond the scope of this paper, but in this section we would like to discuss 
the qualitative features of this deviation. 

To do so, we present an example of a specially useful distribution [123] to probe New 
Physics, namely the invariant mass of the Higgs and vector boson system in Higgs associated 
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production, Myn- as it corresponds to the parton energy, s = My H . In Figure 9 we show 
how the invariant mass distribution changes when dimension-six operators are included. 
Grey and pink distributions correspond to SM and EFT with chw = ~cw = 0.004. The 
width correspond to varying the factorization scale p, = rrih + mz in a range of 2/r and p/2. 
The value of cw is consistent with Runl data and will be explored in the Run2 LHC. The 
value 5 is the deviation respect to the SM central value. 

With the rough estimate given in Eq. 4.2, one can see that for a 2HDM, with A = 3x47r 
the breakdown occurs at around Myu ~ 360 GeV. For the dilaton, in the case of A = 500 
GeV, the breakdown occurs at about 320 GeV. These numbers highlight the importance 
to relate EFT distributions with specific UV models. 

The distributions are the result of a simulation including PDF effects (NN23L01) and 
parton-showering and hadronization performed using an implementation from MCFM [124] 
into P0WHEG [125] at NLO in QCD. We have checked that these results are consistent with 
a parallel implementation in MG5 MC@NL0 [126] using the model implementation in [127] 
into NLO QCD, based on Feynrules [128] and the UFO [129] format. Note that the EFT 
distribution does not grow with y/§, but falls down due to the PDF effects. 

Besides this distribution, one can find others sensitive to New Physics effects, e.g. 
rrijj and A y.jj in Vector Boson Fusion [130] (see also [131]), or Aand pf^ h in di-Higgs 
production [132]. 

5 Discussion and Summary 

During Run2 of the LHC, interpretation of data using the Effective Field Theory (EFT) 
approach will become a standard way to communicate results to theorists, as the translation 
to UV models is more direct than, for example, the use of form factors. 

The approach based on anomalous couplings and the one considered here are related. 
New Physics affects the behaviour of SM particles, inducing anomalous couplings, for exam¬ 
ple involving the Higgs and electroweak bosons. In the context of heavy New Physics, the 
anomalous couplings can be translated into combinations of higher-dimensional operators 
as shown in Tables 5 and 6. 

The use of the EFT approach, however, is limited to heavy new physics, whereas 
current LHC limits strive to reach the TeV region, leading to a delicate balance between 
using as much data as possible and the model-independent EFT approach. The question of 
what LHC data is suitable to constrain EFTs depends both on the UV completion and the 
type of signal one is looking at. Benchmark models allow to address these issues, and also to 
draw correlations with direct searches and other non-LHC sources of indirect constraints 
such as flavour physics and LEP. In this paper we have taken a first step towards this 
program by studying the matching between UV completions encompassing extensions of 
the Higgs sector and the EFT. 

This matching allowed us to study the suitability of LHC data to constrain EFTs. For 
example, one can consider UV completions with large couplings where the validity of the 
EFT is improved, such as the 2HDM in Section 3.2 with large quartic couplings. But even 
for weakly-coupled UV completions, the use of the EFT may be justified as long as the 
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LHC data one uses is restricted to small values of \fi, below the masses of new particles, 
as discussed in Section 4. 

We have shown the matching of the UV theory with the low-energy coefficients in detail. 
This matching is straightforward in the case of the tree-level dilaton exchange (Section 3.3) 
or mixing with a singlet (Section 3.1). For loop-induced dimension-six operators, however, 
the interplay of with higher order terms (dimension-eight operators) needs to be handled 
with care as shown in Section 3.2.2. 

Another advantage of the comparison with UV completions is to address correlations 
among the EFT coefficients which are present in models and reduce the number of free 
parameters in a global analysis. For example, in our benchmarks specific relations among 
operators can be traced back to the (limited) Lorentz structures one can build up from 
scalar fields, see discussion around (3.28). 

Note that we have not discussed the very interesting possibility of CP-violating effects 
from, for example, complex parameters in the 2HDM. This deserves further study, as 
their kinematic distributions and dependence with energy have to be simulated as close as 
possible to the actual cuts applied by the LHC collaborations. 

Note added: As this work was being submitted for publication, we became aware of [136], 
which also discusses UV completions of Higgs Effective Field Theory via extended Higgs 
sectors finding results similar to ours. We stress that both works are complementary: we 
focus on the connection between the EFT and UV models as a way to assess the range 
of validity of the EFT, whereas [136] discusses the connection between the EFT and UV 
models in the context of a systematic framework for the obtention of the EFT Wilson 
coefficients from an arbitrary UV theory. 
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A 2HDM Results and Conventions 
A.l Scalar Potential and Mass Spectrum 

Let us recall the scalar potential for a 2HDM with a softly-broken Z 2 -symmetry in the CP- 
conserving case, given by (3.25). After EW symmetry breaking, the scalar mass eigenstates 
can be written in terms of the original fields in (3.26) as 

G ± = cos (3 ipf + sin j3 H ± = — sin f3 tpf + cos (3 ipf 

G° = cos (3 7/i + sin f3 r )2 A 0 = — sin (3 rji + cos (3 772 (A.l) 

h = — sin a h\ + cos a /12 H° = — cos a h\ — sin a /12 
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with i? , A 0 , H°, h being the physical states of the theory and G , G° the Goldstone 
bosons from the breaking of EW symmetry. Requiring that the scalar potential be bounded 
from below yields 


Ai > 0 , A2 > 0 , A3 > —a/A1A2 , A3 + A4 — | A5I > — \J A1A2 • (A.2) 

Obtaining the correct EW vacuum as the minimum of the scalar potential (3.25) imposes 
the relations 


= A1C2 + A34B-J), 

h 2 2 = m 2 4- v ^(\ 2 s}+\ 3 , 5 4), 


(A.3) 


with A 345 = A 3 + A 4 + A 5 , sp = sin/3, cp = cos (3 and M 2 = /i 2 / spcp. Similarly, the quartic 
couplings Ai _5 in (3.25) may be re-expressed in terms of the masses of the physical states 
m H ±, m^o, m H 0 , m-h, the mixing angles a , f3 and M 2 as 


Ai = 


A 2 — 


v 2 c 2 ^ 


V 2 S 2 p 


(-M 2 s 2 p + m\sl + m 2 H0 c 2 ) , 
(—M 2 c 2 + m^o 4) , 


A 3 — —5 


— M 2 + 2m 2 H ± + (m^o - m 2 ) 


A 4 = ^ (M 2 + <0 ~ 2m#±), 

A s = 4 (Af 2 - m^o) • 


(A.4) 

(A.5) 

(A. 6 ) 

(A.7) 

(A. 8 ) 


We stress that since the masses of the physical states m H ±, m^o, m H 0 , m?, and the mixing 
angles a, /3 are obtained upon EW symmetry breaking, the relations (A.4-A.8) only hold 
in the EW broken theory. 

We now perform an 517(2) rotation from the field basis 4>i, 4>2 to a field basis Hi, H 2 
in which only H\ takes a vev: (H\) = , (H 2 ) = 0. This rotation is given by 


H 1 = cp 4>i + sp $2 H 2 = — sp 4>i + cp 4*2 


As seen from (A.l), for a = f3 — 7 r/ 2 , we obtain after the field rotation 


Hi = 4>sm 


o* \ 

v-\-h-\-iG° I 1 

V2 ) 


( H+ \ 

H 2 — | H°+iA° I > 

V y/2 ) 


(A.9) 


(A.10) 


This is the alignment limit cp- a = 0, in which H\ corresponds exactly to the SM Higgs 
doublet. Away from this limit there is mixing of the neutral CP-even physical states h, H° 
between Hip. 
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The scalar potential for H\ 2 is given by 


Vt ree (Hi, H 2 ) = ft \Hft + ft \H 2 \ 2 - ft h\h 2 + H.c. 
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T Ag 
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\H 2 \ 2 h\h 2 + H.c. 


(A.ll) 


The modified mass parameters ft, ft, ft and quartic couplings Ai_s are expressed in terms 
of m 2 H ± , m 2 A o, m 2 H o , rft, M 2 , a and (3 as 


~2 1 r 2 2 1 22 1 , n 

[ m H a c p_ a + rn h Sp_ a \ < 0 


ft = M 2 - - (m 2 H 0 - ml) ^ + m 2 H0 s}_ a + m 2 h c}_ c 
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(A.12) 
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(A.13) 


The relations ft = —2 \\ v 2 and ft = v 2 are necessary to obtain correct vev for H[ 2 
upon minimization of (A.ll). The relations among the masses of the new scalar states may 
be written in terms of xq = m 2 H0 / m 2 H± and xa = m 2 A o / r fti± as 
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In the alignment limit, the above expressions (A.12-A.13) simplify considerably, yielding 
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(A.17) 


The interaction vertices ghH+H~i 9hA°A° and ghH°H° are in this limit proportional respec¬ 
tively to A 3 , A 3 + A 4 — A 5 and A 3 + A 4 + A 5 . 


A.2 Electroweak Precision Observables 


We turn now to the discussion of EW precision constraints in the 2HDM. The non-SM 
contributions to the oblique parameters S and T are given by (3.50-3.51) (the general 
expressions for the various oblique parameters in models with an arbitrary number of 
scalar doublets may be found in [113]), with the functions Fa,Bi Ga,b,Ci Ga,b being 
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(A-19) 


Finally, we stress that the potential (A. 11) preserves custodial symmetry in the limit A 4 = 
A 5 (see e.g. [133]), which from (A. 14) corresponds to m^o = m H ±. To see this, instead of 
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expressing the 2HDM scalar potential in terms of Hi 2 , we can introduce the 2x2 matrices 
= (iooHl , H\), $2 = # 2 )- The scalar potential for the 2HDM then reads 
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(A.20) 


Both $i j2 transform as bi-doublets of a global symmetry SU{2)l x SU(2)r: <f>, —>• L<& t R 
(with L E SU(2)l and f? E SU(2)r). The potential (A.20) is then invariant under a 
custodial SU(2 )l x SU(2)r in the absence of the term proportional to A 4 — A 5 . 


B Effective Theory Matching in the EW Broken Phase 


We now relate the Wilson coefficients from the D = 6 operators for the SM SU(2)l xU(1)y 
invariant effective field theory given in (2.2), (2.5) and (2.6) to the SM effective Lagrangian 
after EW symmetry breaking. The complete set of relations from the full D < 6 effective 
Lagrangian (2.1) may be found in [26]. In the following we present the generic relations for 
the 3- and 4-point interactions involving Higgs/gauge bosons. The CP-conserving 3-point 
interactions involving at least one light Higgs scalar h are obtained from Tsilh and read 
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- ^ Z^F^h - g^l Z v d^h . 


(B.l) 


where we have introduced abelian field-strength tensors W )iU . Z^ v and F^ for the W-boson, 
Z-boson and photon respectively, and is still the non-abelian gluon field-strength 
tensor. The CP-conserving 4-point interactions involving at least one light Higgs scalar h 
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similarly read 
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We note that the 4-point interactions involving one Higgs h and three gluons (relevant, 
e.g. for Higgs production in association with a jet via gluon fusion) are already taken into 
account by the g hgg term in (B.l). In a similar fashion to (B.l) and (B.2), the CP-odd 3- 
and 4-point interactions involving at least one light Higgs scalar h are obtained from Cqv, 
yielding 
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(B.3) 


Finally, the 3- and 4-point gauge boson self-interactions receive CP-even contributions from 
both £silh and jC-g (they also receive CP-odd contributions from £cp> which we do not 
include here). We adopt here the parametrization of triple gauge couplings (TGCs) from 
[10, 134], which yields 
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with g[ = gf = Kry = n z = 1 and A 7 = \ z = 0 for the SM. U (1)em gauge invariance imposes 
g\ = 1 , while the (spontaneously broken) SU(2)l x U(1)y gauge invariance combined with 
the restriction to D < 6 effective operators in £Eff lead to the relations [10] 
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C^h Couplings vs SU(2 )l x U(1)y (D < 6) Wilson Coefficients 
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Table 5. Relations between the different couplings appearing in £ 3 ^, £ 4 ^, and the Wilson 
coefficients for the D < 6 effective operators in £es- 


These degeneracies imply that there can be at most three independent contributions 
to anomalous TCGs from £silh and Cq■ The operators proportional to gf, k z and « 7 in 
(B.4) contain one derivative and are obtained from £silh> which then gives rise at most 
to two independent contributions to anomalous TGCs. The last two operators in (B.4) 
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C 3V Couplings vs SU(2 )l x U(1)y (D < 6 ) Wilson Coefficients 
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Table 6. The Trilinear and Quartic Gauge Couplings appearing in £ 3 ^ and £ 4 ^ as a function of 
the Wilson coefficients for the D < 6 effective operators in £es- 


contain three derivatives and are obtained from £< 7 , which then gives rise to a sole linearly 
independent contribution to anomalous TGCs. This is clear for the operator c 3W 0^w , as it 
is constructed from three gauge field-strengths and only their abelian parts can contribute 
to TGCs (after EW symmetry breaking it directly leads to the last two operators in (B.4) 
with A 7 = A z ). The operator c 2 wOdw also leads after EW symmetry breaking to those 
operators with A 7 = A z , but also produces operators of the form 


ge ijk {UW^WlW* 


(B. 6 ) 


These operators however either vanish or reduce to operators from (B.4) for on-shell gauge 
bosons (as well as for the virtual photon). Therefore, the parametrization of TGCs (B.4) 
holds. Moreover, we can safely disregard the contribution of c 2W to A 7 = A z in the following 
discussion, since this operator can be re-casted as (very constrained) fermionic currents and 
the operators involving the Higgs we have already considered [135]. For the case of quartic 
gauge couplings (QGCs) we have 
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where again = fff = 9]. Z = 1 an d -W = A 7 w = A 7Z = Aiy Z = 0 corresponds to the 

SM case. 


The correspondence between the Wilson coefficients Cj for the effective operators from 
(2.2), (2.5), (2.6) and the various couplings arising after EW symmetry breaking in (B.l- 
B.4), (B.7) is shown in Tables 5 and 6. In the next section, we compute the values of these 
couplings by direct matching in the broken EW theory, for the specific setup of Benchmark 
A in 3.2.1: the 2HDM in the alignment limit cp_ a = 0. 


B.l A Specific Example: 2HDM in the Alignment Limit cg_ Q = 0 

Here we will obtain the values of the various couplings from (B.l) for the 2HDM in the 
alignment limit cp- a = 0, by direct matching in the broken EW theory. We focus on the 
Vfj,V u h interaction vertices, which receive contributions from loops of scalars H^,A 0 ,H 0 
(these are shown in Figure 10 for the case V) ( = Z^). We note that loops involving only 
SM particles cancel in the matching between full and effective theories and may then be 
disregarded, while loops involving both SM particles and H±, A 0 , H° scalars vanish in the 
alignment limit. 


h 


h 



Z u 




Figure 10. 1-loop Feynman diagrams involving the H ± ,A°,H° scalars and contributing to the 
Z^Z V h vertex in the alignment limit cp- a = 0 


We compute the 1-loop contributions to V^V u h and expand the result in powers of the 
4-momenta of the Z-bosons pi,p 2 - By truncating the expansion at order 0(p 2 ), we obtain 

AC = A C\y + A Cz + A£ 7 (B.8) 
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with the last three terms in both (B.9) and (B.10) corresponding to e.o.m. -vanishing terms, 
which are nevertheless generated by the off-shell 1-loop corrections. Performing an expan¬ 
sion to linear order in 1 — m 2 H0 /m 2 H± = 1 — xo and 1 — m 2 A0 /m 2 H± = 1 — xa yields for 
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with go = ghH°H°/ v i 9A = 9hA°A°/ v an( 4 9+ = 9hH+H~ t v i and the trilinear scalar couplings 
9hH+H ~> 9hH°H° and 9hA°A° given in the EW broken theory with cp- a = 0 by 
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By means of (B.15), the relations (B.12 - B.14) may be written as 
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Performing a similar expansion in A Cz then yields 
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while for A£ 7 we obtain 
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Use of (B.15) in (B.19 - B.22) yields 
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Altogether, these results show that g^yy ^ ~ v 4 /\x\ at least, and only receive contributions 
from D = 8 effective operators. For the D = 6 effective operators, we use the leading order 
in (B.16 - B.18) and (B.23 - B.27), and through the relations from Table 5 we get 
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which reproduce the results from (3.28) upon the substitution m 2 H ± —> jj 
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